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Electromagnetic resonant cavities form the basis of many tests of Lorentz invariance involving 
photons. The effects of some forms of Lorentz violation scale with cavity size. We investigate 
possible signals of violations in the naturally occurring resonances formed in the Earth-ionosphere 
cavity. Comparison with observed resonances places the first terrestrial constraints on coefficients 
associated with dimension-three Lorentz-violating operators at the level of 10"^" GeV. 
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I. INTRODUCTION 

Modern versions [1] of the classic Michelson-Morley [2] 
and Kennedy-Thorndike [3] experiments are among the 
most sensitive tests of Lorentz invariance, the symme- 
try behind special relativity. Typically these experiments 
search for minute changes in the resonant frequencies of 
electromagnetic cavities. High quality factors allow for 
precise tracking of the frequencies, giving extreme sensi- 
tivities to possible deviations from perfect Lorentz sym- 
metry. However, the effects of some forms of Lorentz 
violation increase with wavelength. As a result, high sen- 
sitivities may be achieved using very-low-frequency reso- 
nances, such as those that naturally occur in the Earth- 
ionosphere cavity, despite their relatively low quality fac- 
tors. 

In this work, we consider signals of Lorentz violation 
that would appear in Schumann resonances, the lowest- 
frequency standing waves that form in the atmosphere 
[4, 5]. We obtain conservative bounds by comparison 
with observations [6]. The Earth's surface and ionosphere 
form a cavity of immense size, leading to resonances with 
very long wavelengths. The lowest-frequency resonances 
have wavelengths that are comparable to the circumfer- 
ence of the Earth and have frequencies as low as 8 Hz. 

Violations of Lorentz invariance are described by the 
Standard-Model Extension (SME) [7, 8]. The SME is 
a theoretical framework that provides a basis for many 
experimental and theoretical studies of Lorentz viola- 
tion [9, 10], including those involving atoms [11, 12], 
hadrons [13, 14], fermions [15-18], the Higgs boson [19], 
gravity [20], and photons [1, 21-25]. In addition to 
resonant-cavity experiments, searches for Lorentz viola- 
tion in photons using the SME approach include astro- 
physical searches for vacuum birefringence [22-25] and 
dispersion [24, 25]. The goal of the SME is the charac- 
terization of all violations of Lorentz symmetry that are 
consistent with known physics using effective field the- 
ory [26]. While motivated in part by the possibility of 
spontaneous symmetry breaking in strings [27, 28] , it en- 
compassed violations with other origins [29-36]. Much 
of the work on Lorentz violation has focused on the min- 
imal SME, which includes operators of rcnormalizable 
dimension in a flat spacetime. While nonrcnormalizable 
operators and curved spacetimes are of general interest 
[8, 20, 23, 24, 30, 37], Schumann resonances are particu- 



larly sensitive to the dimension-three CPT-odd Lorentz- 
violating operators of the minimal-SME photon sector. 

In cavities, Lorentz violation can introduce frequen- 
cies that depend on the orientation of the cavity, sig- 
nalling rotation violations, and dependence on velocity 
resulting from boost violations [25]. The quantity de- 
termining their sensitivity is the dimensionless fractional 
frequency shift Svjv. To date, cavity experiments have 
focused primarily on one class of violations, namely the 
dimension-four C PT-even Lorentz-violating operators of 
the minimal SME. The coefficients associated with these 
operators (/cj^^)"'^''''^ are dimensionless. Therefore, di- 
mensional analysis suggests frequency shifts of the form 
~ (kp^)"^'^^, which implies little or no depen- 
dence on frequency. Consequently, there is little advan- 
tage to using low-frequency resonances. In contrast, the 
coefficients associated with the dimension-three opera- 
tors (fc^^)^ have mass-dimension one. Therefore, we 
naively expect shifts in frequency that depend on the 
ratio {k%)^lu. Given that u ~ 10-^3 GeV for Schu- 
mann resonances, we naively expect sensitivities on the 
order of 10"^^ GeV to fc^]^ coeflticients, assuming at least 
order-one sensitivity to 8vjv. While not as sensitive as 
birefringence tests [24], the bounds obtained here repre- 
sent the first terrestrial bounds on the dimension-three 
operators, providing a valuable check on existing astro- 
physical constraints. 

The structural outline of this paper is as follows. Sec- 
tion II provides some basic theory behind our calcula- 
tion. In Sec. Ill we derive modified wave equations and 
describe a numerical method of determining the effects 
of Lorentz violations on Schumann resonances. The re- 
sults of our calculation are discussed in Sec. IV. Unless 
otherwise stated, we use the notation and conventions of 
Refs. [7, 24]. 



II. BACKGROUND 

In this section we discuss the theory behind the calcu- 
lation of the Earth-ionosphere resonances in the presence 
of dimension-three Lorentz violations. We begin by dis- 
cussing the conductivity of the atmosphere and the Schu- 
mann resonances in the usual case. We then review the 
modified electrodynamics including the CPT-odd oper- 
ators of the minimal SME. 
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A. Conductivity profile 

Resonances in the Earth- ionosphere cavity are excited 
by a number of man-made and natural phenomena, light- 
ning being a primary source. The surface of the Earth 
forms the lower boundary and, in our calculation, is 
treated as a perfect conductor. The ionosphere forms 
a lossy upper boundary with a finite conductivity pro- 
file that increases with altitude. The conductivity of 
the lower atmosphere can be approximated by a "knee"- 
model that separates into two layers with exponentially 
increasing conductivity [38]: 

{00 r < R , 

aa cxp ^ R<r <ro , (1) 
(70 exp To <r , 

where r is the distance for the center of the Earth, and 
R ~ 6400 km ~ 3.2 x 10^2 GeV"^ is the Earth radius. 
The lower layer is dominated by positive and negative 
ions. The upper layer approximates the bottom of the 
ionosphere, which is dominated by free electrons. The 
knee radius vq is the boundary between the layers, and 
(To is the conductivity at this transition. In units where 
c = ?l = Co = (Uo = 1) we adopt values ro = 1.009-R, 

= O.OOTi?, 4 = 0.0005i?,, and ctq = l.Si?-^ in the 
numerical calculations that follow. This profile yields 
frequencies and quality factors that closely match the 
observed resonances in the Lorentz-invariant limit. 

In the conventional case, both transverse-magnetic 
(TM) and transverse-electric (TE) modes may be excited 
in the cavity. However, TE modes oscillate in the radial 
direction, implying wavelengths comparable to the lieiglit 
of the ionosphere, yielding frequencies in the kHz range. 
In contrast, the lowest-frequency TM modes vary little in 
the radial direction but form standing waves that encircle 
the Earth. As a result, the wavelengths are comparable 
to the Earth's circumference, yielding frequencies as low 
as 8 Hz. Our goal is to understand the effects of Lorentz 
violations on these low-frequency Schumann resonances. 

B. Modified electrodynamics 

The lagrangian governing electromagnetic waves, in- 
cluding dimension-three Lorentz-violating operators, is 
given by [7] 

£ = -iF^.F^-^ + ie'^^^''(fciF)«^AF^. , (2) 

where Ffj_„ = d^Ai, — d^Afj, is the field strength. The re- 
sulting theory preserves the usual gauge symmetry, but 
violates CPT invariance. Lorentz and CPT violations 
arc controlled by the constant coefficients {k^'^p)^, which 

include a pseudoscalar (fc^]^)o and pseudovector fe^^. 

The (fc^^)^ coefficients are assumed to be constant, 
which leads to energy-momentum conservation. More 
generally one can consider Lorentz-violating backgrounds 



with spacetime variations. These types of violations are 
particularly important when considering Lorentz viola- 
tions in curved spacetimes [8, 20]. However, in scenarios 
where violations originate in the moments shortly after 
the big bang, it is likely that any variation has expanded, 
leading to little fluctuation over experimentally relevant 
time and length scales. Consequently, the idea presented 
here probes {k\p)K, in our local neighborhood. Astro- 
physical tests rely on light that has propagated across 
much of the visible universe. So they test Lorentz in- 
variance over much larger scales and could be drastically 
affected by variations in {k\p)n- Therefore, while current 
bounds from astrophysical searches for birefringence cur- 
rently lie at the 10"^^ GeV level [24], terrestrial tests 
provide an important complementary set of constraints. 

The equations of motion resulting from Eq. (2) pro- 
vide Lorentz-violating inhomogeneous Maxwell equa- 
tions. We arc primarily interested in harmonic solutions 
and consider electric and magnetic fields of the form 
E{t) = E{uj)e-^\B{t) = B{uj)e-^'^K Together with 
the usual homogeneous equations, we arrive at a Lorentz- 
violating electrodynamics with a modified Ampere law 
and a conventional Faraday law: 

V xB + iojE - 2 {k'-^l)o B + 2 k^^l x E -aE = , 

V X E- iujB = . (3) 

Here we have included the usual source current and con- 
ductivity terms. These conventional source terms result 
if we assume usual coupling to matter. The source cur- 
rent is not needed in determining the resonances but 
could be used in modeling the effects of individual light- 
ning strikes. The conductivity term is necessary to ob- 
tain realistic frequencies and quality factors, which are 
affected by the profile of the upper boundary and the 
losses it introduces. 



III. CALCULATION 

Our goal is to calculate the resonant frequencies that 
result from the modified Maxwell equations. We begin by 
expanding into spherical harmonics and deriving modi- 
fied wave equations. A numerical calculation is used to 
estimate the resulting resonant frequencies. The results 
are discussed in the following section. 



A. Wave equations 

We begin our search for resonances by first expressing 
the modified Maxwell equations (3) in spherical coordi- 
nates using the helicity-basis and identities discussed in 
Appendix A. This involves writing the Maxwell equa- 
tions in covariant form, then using relation (A5) to ex- 
press them in terms of d/dr and covariant angular- 
momentum ladder operators J±. Dropping the source 
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current, the result is 



= J+B_ + J_B+ + {uj + ia)rEr + 2ir{k%)o^r 

= ±§^rB± - J±Br + (w + ia)rE± + 2ir{kfp)oB± 

± 2r{kfp)rE± T 2r{kfp)±Er , 
= J+E_ + J_E+ - ujrBr , 



u = J+E- + J-E+ — (jjrBr , 
= ±§^rE± - J±Er - wrB± , 



(5) 
(6) 
(7) 



where Er, Br are the radial field components, and E±, 
B± are negative/positive holicity components, as dis- 

(3) 

cussed in the appendix. The components of k\p are 

defined by (fc^^)a = e.a ■ f^AF^ where arc the hchcity 
basis vectors. The components E± and B± arc referred to 
as spin- weighted functions with spin-weight- (±1), while 
Er, Br have a weight of zero. Spin weight and helicity 
are equivalent up to a sign. Eqs. (4) and (6) provide two 
scalar relations, while (5) and (7) have a spin weight of 
±1. 

We can expand the field components in spin-weighted 
spherical harmonics. These provide the generalization of 
the familiar spherical harmonics to spin-weighted func- 
tions. The expansion takes the form 



(8) 



E± 



E 



2 r 



(te^''^^ 



iE^l^^)±rYjm, (9) 



j(j-H) 1 
2 r 



{TB. 



(IB) 



(10) 

jm ) (11) 



where E] 



(OE) pilE) „(0B) „„j d(1E) 

'jm ' jm ' jm ' jm ' -^jm ' "^^^^ jm 

are r-dependent field coefficients. They are associated 
with total-angular-momentum eigenmodes and have E- 
type parity, (— l)-*, or B-type parity, (—1)-'+^. The 
+ l)/2 and 1/r factors in the expansions are for 
convenience. 

Using these expansions, we can express the Maxwell 



equations in terms of eI°^\ E^l^\ EfJ^>, B^^ 



(OS) 



B 



(IB) 



jm ' jm ' 

and Bjl^\ First, using the ladder operators Jip to 

lower/raise the spin-(±l) relations (5) and (7) wc ar- 
rive at six scalar Maxwell equations. We then use the 
spherical-harmonic expansions of the fields and the or- 
thogonality relation (A3) to get relations between the 
six expansion coefficients. Some algebra yields three im- 



parity equations and three S-parity equations: 

= j{j + l)B^l^^ - iiuj + ia)rEf^^ + rJCf^^ , (12) 



" dr^jm 



.i[u^ + ia)Ef^^ + ^l6^\ (13) 



(4) = ri,<)-<)+^(c. + ^a)r<) 



I r ir(iB) 
lU+T)'^jm ' 



= JiJ + l)E^lf^+^u;rBly 

" dr-'^jm ^ ''^-^jm ' 



^ — ' dr ^jm 



E 



(OB) 



lUjrB 



(IE) 



(14) 

(15) 
(16) 
(17) 



where the Lorentz- and CPT- violating contributions 
have been collected into the field combinations 



2\k 



(3) 



Am 



AfKJ + '^)'^{j+l)mE\^+i)m 



(18) 



(IB) 



2{kfl)o]{j + l)B 



(IB) 



,(3) I 
,(3) 



p(lS) 
^jm 



,(3) , 
"AFl 



■^jm 



2j - 1)C,,^eI]'^1 



(i-l)m 

(IB) 



+ 2\k^p\{j - 2)C(j_|_i)mi?^^.^-^^ 



(19) 



icf^^ = -2(fc(f) )oj(j + 1)5]™ ' - 2i\k^2UmE; 
-2|fc(f],|(j + l)C,™£;;°fj)„ 



(IB) 
jm 



+ 2|fe(f),|iC(,+i)„£;[°f))„ 
+ 2|fe(f],|(/ + 2j-l)C,„,i;;jfj 
+ 2|fe(f],|(/-2)C(,.+i)„£g2), 



(20) 



where Cjm = -^Z {j^ — ■m?)/{Ap — 1). Here we take the 
angular-momentum quantization axis along the direction 

oik\p. Note that Eqs. (12)-(14) correspond to the mod- 
ified spherical Ampere law and Eqs. (15)-(17) are the 
usual Faraday law. 

In the conventional case, where all coefficients for 
Lorentz violation are zero, rotational symmetry implies 
that resonances are eigenmodes of angular momentum 
with definite j and m values. The symmetry also im- 
plies degeneracy in m. So the different resonant frequen- 
cies correspond to different values of the total-angular- 
momentum index j. Setting A^]^'* = ^j>,f ■* = ^jLf = 
in Eqs. (12)-(17), we also note that the Lorentz-invariant 
case splits according to parity. The J5-parity resonances 
correspond to the high-frequency TE modes, while low- 
frequency TM Schumann resonances are the E'-parity 
modes. 

Allowing for Lorentz violations, the new symmetries of 
the system lead to several generic predictions. Rotational 
symmetry is preserved in the event that we have only 
isotropic violations associated with coefficient (fcj^p.)o- 
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This implies that the indexing and degeneracies of the 
modes is the same, but the frequencies may change. In 
contrast, the vector k\p breaks the usual degeneracy. 
The system remains symmetric under rotations about 
^AF- expect resonances that are eigenmodes of 

these rotations with eigenvalues m, as usual. However, 
those coefficients break spherical symmetry, implying the 
index j is no longer associated with resonances. As a re- 
sult, the usual 2j + 1 degeneracies should break, yielding 
modes with definite m but indefinite j. Consequently, 
we expect two types of effects that would signify pos- 
sible Lorentz violation. One is a split of degeneracies 
leading to additional resonant frequencies. This results 



from anisotropic violations. The other effect is a shift in 
frequencies that may result from either anisotropic and 
isotropic violations. 

The above first-order differential equations can be re- 
duced to second-order modified wave equations. We be- 
gin by using the Faraday law, Eqs. (15)-(17), to elimi- 
nate the magnetic field. We also use Eq. (12) to elimi- 
nate the electric field component Ej^^ in favor of Ejl^^ , 

E^l^\ and icf^\ The result of this process is three cou- 
pled equations relating the three sets of field components 

Ejm ' ^jm ' and Ao^„ . 



n_9^ I P^i (d iv+ic7 \a pirn „2p{m ip^i 1 ir(Og) I o„-|ij3) I ma- ^(0^) 

0//,(3)y a pjiB) _9|r.(3) I P"" pirn o|. (3) l m d pjlE) 

^V^-^F/O dr-^jm ^I'^'AfI {u;+ia)j{j+l) ^jm ^ ^I'^AfI {u;+ia)r dr^jm 

T-^'I^-AfI {u+icr)j(3 + l) -^0-1)™ ^'I'^'AfI (w+io-)j(j+l) ^(j+l)m^ \^^) 

n-d^ 2p{lB) 017.(3) I rnuj p(lB) (3) y {u+ta)uj 3 p(^E) r,- ( A3) cj ic(OE) 
^—W^^jm +Pj-^jm ^ '^\'^AF\j{j + l)^jm +^l'^AFiO pi W-^jm + ^H'^AfJo 

o„-|/-.(3) | (j"-H2j-l)C,„c^ pjlE) o-|ij3) | (/-2)Co + i)„c^ (IF) o -I |J3) , (j-l)C,„.t^ a ^(IF) 

- ^»I«AfI JU+T) ^(i-l)m ^^I'^AfI JJj+T) ^ij+l)m ' ^^I^AfI pj_^r 9F^{j-l)m 

I 0|l.(3) I CjmO'^ 

+ ^l«AFl7^^(j-l)m + ^*I«AfI aF-^(j-M)m '^l«AFl(j + l)p^^^^%+l)m ' l^^j 

- 2^(e)o^<^ - 2i\k^lUmEll^' 

- 2\kfp\{j - 1)Cj^e'^%^ + 2\k^2l\{j + 2)C(,.+i)„i?;;^])„ , (23) 



= 



where we define p| — uj{uj + ia) — j{j + 1)/^^. Note 
that we could use Eq. (23) to ehminate icf^^. How- 
ever, for simplicity, we treat /Cj^"* as a dynamical field 
on equal footing with Ejl^^ and Ej^^ . The field compo- 
nents Ejl^^ and Ejl^^ correspond to the transverse part 
of the electric field, which vanishes at the surfaces of a 
perfect conductor. This implies /cj^'' vanishes on the 
surfaces as well. So we take e'^1^\ -E^if' , and IC[^^ as 

J III' ' J III' J III' 

independent fields that vanish at the boundaries of the 
cavity. 

B. Numerical frequencies 

We calculate the resonant frequencies that result from 

the modified electrodynamics by considering discrete 
radii = R + Sr{n+ |), where n is an integer. Defining 

discrete field coefficients at these points, E^^j^, E^^j^, 

and )Cn^' a^*^ using discrete derivatives, wave equa- 
tions (21)-(23) can be written in the form of an infinite- 
dimensional matrix equation. Discrete resonances cor- 



respond to frequencies where nontrivial field configura- 
tions exist. These can be estimated by truncating the 
matrix at finite index values and searching for uj where 
the truncated matrix is singular. These u) are complex 
in general. The real parts give the resonant frequencies 
V = Rea;/(27r), while the ratios of the real and imaginary 
parts determine the quality factors Q = — Rea;/Ima;/2 
of the modes. 

In the event that ik^^l)^, = 0, Eqs. (21)-(23) reduce 
to two wave equations, one for B modes and one for the 
Schumann E modes. However, nonzero (fc^^)^ coeffi- 
cients mix the two parities, and resonances will no longer 
possess definite E ov B parity. We also note that no mix- 
ing of fields with different m values occurs, but mixing 
across j values results from the vector k)^p, as expected. 
As a result, all resonances have definite m values, and m 
may be fixed in the calculation. 

To determine the resonances, we take 100 different 
r values, n — 0, . . . 99, uniformly spaced between 7? 
and l.Oli?. This corresponds to a penetration of about 
O.OOlii ~ 6.4 km into the highly conductive ionosphere. 

For fixed m, we create a matrix inchiding terms corre- 
sponding to these n values and the ten lowest j values 
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FIG. 1: Calculated resonant frequencies and quality factors 

for three different values of [(fe^^)o|- Circles represent the res- 
onances. The Lorentz-invariant limit is shown with x symbols 
for comparison. The shaded regions indicate the correspond- 
ing widths at half maximum for the three lowest resonances 
in the Lorentz-invariant case. 



that are relevant, j > \m\. The result is a square matrix 
with dimension 100 x 10 x 3. We use a row-reduction 
method to determine its determinant for different values 
of u) and search for roots. 



IV. RESULTS 

While any combination of coefficients for Lorentz vi- 
olation is possible, for simplicity, we next consider the 

(3) (S^ 

effects of (fc^^)o and Ife^^^'l separately. We first consider 

a nonzero {k)^p)o coefficient. Figure 1 shows the three 

lowest-frequency resonances for three values of |(fc^^)o|. 
The three resonances shown correspond to j = 1,2,3 
and are degenerate in m. The results are independent 

/ON 

of the sign of {k\p)o due to the symmetry of this case. 
The figure also shows the calculated resonant frequency 
and Q factor for the Lorentz-invariant limit. These 
are in good agreement with the observed resonances of 
(i/,Q) = (7.8 Hz, 4.0), (14.1 Hz, 4.5), (20.3 Hz, 5.0) [6]. 

From the figure we see that values of {k\p)o on the 
order of 10-^^ GeV significantly affect both the resonant 
frequencies and Q factors. In particular, values of 4 x 
10-^1 GeV drastically alter all three modes. We therefore 
adopt a conservative limit of \{k^^p)T\ < 4 x 10-^^ on the 
time-like part of (fc^^)^ hi the standard Sun-center frame 
described in Ref. [25] . This translates to a bound of 



|fc(y)ool < 14 X 10-21 
on the spherical coefficient from Ref. [23]. 



(24) 



FIG. 2: Resonant frequencies and Q factors for three values 
of The |fe^^| — case (x symbols) is shown for 

comparison. The shaded regions indicate the corresponding 
widths at half maximum for the three lowest resonances in 
the Lorentz-invariant case. Each plot includes resonances for 
index values m = (stars), m = ±1 (filled/empty circles), 
m = ±2 (filled/empty squares), and m = ±3 (filled/empty 
triangles). 



Note that the above bound is about two orders of mag- 
nitude larger than the naive prediction. This can be un- 

derstood from the fact that the {k\p)K coefficients mix 
E- and iJ-parity modes. The i3-parity resonances have 
frequencies that are at much higher frequencies. This 
leads to a seesaw effect in the Maxwell equations that 
suppresses perturbations in the resonances. As a result, 
relatively large mixing in the wave equations must be 
present for significant changes to manifest in the low- 
frequency modes. 

Small changes to the conductivity profile can lead to 
large changes in the resonances, implying that our confi- 

dence in the bound on {kj^p)T is somewhat weakened by 
our knowledge of the atmospheric conductivity, which is 
a complicated and variable system. Much cleaner bounds 

(3) 

can be placed on the pseudovector part k^^p since it 
leads to a breakdown of the usual 2j + 1 degeneracy 
among modes with identical j eigenvalues. These bounds 
too are complicated by imperfections in our conductivity 
model, including the missing day-side/night-side and po- 
lar asymmetries that are present in the real atmosphere 
[5]. These conventional anisotropics can also break the 
degeneracies, but they would only add to the effects we 
are bounding. We therefore neglect them in our analy- 
sis. However, these features could be significant in more 
detailed studies involving field configurations. For exam- 
ple, local fields may experience variability that includes 
daily and annual fluctuations from conventional physics. 
Extracting a sidereal dependence caused by the rotation 
of the Earth with respect to the fixed k)^p vector might 
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be possible but is beyond the scope of this work. 

Figure 2 shows the resonances for three different val- 
ues of We notice shifts in the frequencies and Q 
factors as well as the expected 2j + 1 splitting of the res- 
onances. In particular, values of = 8xl0~^^ GeV 
lead to new resonances separated by frequency intervals 
comparable to the resonance widths. These multiple res- 
onances would be evident in the data if they existed. 
Therefore, we take \k^2l\ < 8 x IQ-^i GeV as a conser- 
vative bound. The relation to the spherical coefficients 
is given by \kf^\ = (elfc}^) ^ + gj^W ^^|2) [34]. 

So our bound leads to two constraints on spherical coef- 
ficients for Lorentz violation of 

\kf^)^^\ < 12 X 10-21 GeV, 

|fc(^nol < 16 X 10-21 GeV. (25) 

These completely bound the vector (j = 1) dimension- 
three Lorentz-violating operators. Again, these con- 
straints are less stringent than the naive estimate. 

V. DISCUSSION 

In this work, we used Schumann resonances in the 
Earth-ionosphere cavity to place bounds on the order of 
IQ-^o GeV on CPT-odd j = 1 coefficients of the min- 
imal SME. Similar bounds are places on j = scalar 
coefficients assuming the actual conductivity profile of 
the atmosphere is not significantly different from our 
model profile. These bounds constitute the first terres- 
trial bounds on dimension-three Lorentz-violating oper- 
ators. 

While not as sensitive as astrophysical searches for 
vacuum birefringence, the techniques used in this work 
test Lorentz invariance in our local neighborhood, giving 
a bound on coefficients over solar-system length scales. 
In contrast, astrophysical tests probe Lorentz violation 
over cosmological scales and may be obfuscated by space- 
time variations or domains in the Lorentz-violating back- 
grounds. So local tests play an important role in our 
search for new physics. 

Laboratory-based experiments may also provide local 
tests of Lorentz invariance. Current cavity experiments 
utilize high Q factors that allow for sensitivities to 5v/v 
on the order of parts in 10^^ or better [1]. This suggests 
improved bounds may be possible in laboratory experi- 
mcnts. A rough estimate yields sensitives to {kj^^p)^ from 
about 10-2^ GeV in optical cavities to around 10-'^° GeV 
in lower- frequency microwave cavities. 

Future studies involving Schumann resonances may be 
able to improve on the above bounds. Precise tracking of 
the resonances may allow for sidereal searches that would 
indicate rotation violations from k)^p. Also, boost vio- 
lations from (fc^^)^ coefficients can lead to annual de- 
pendences that may be discernible. Regardless, the con- 



straints obtained here demonstrate the potential of res- 
onator experiments as tests of dimension-three Lorentz 
violations in the atmosphere and in the laboratory. 



APPENDIX A: SPIN WEIGHT 

The problem addressed in this work is most naturally 
solved in spherical coordinates. Here we use a helicity- 
based system and covariant-angular-momentum opera- 
tors. In this appendix, we summarize some the key iden- 
tities used in the calculation of the modified wave equa- 
tions. A fuller discussion of these methods will appear 
elsewhere [39]. The technique is based on a decompo- 
sition into total angular momentum J and helicity. A 
type of tensor spherical harmonics called spin-weighted 
harmonics gYjm [40] provide orthonormal sets of eigen- 
functions of these operators. The index s labels the spin 
weight, which up to a sign is equivalent to helicity. 

The method starts by defining helicity basis vectors, 
e± = = -^{ee ± ie^), e-r = e*", where = x/\x\ is 
the radial unit vector, and €.$ and e<^ are unit vectors as- 
sociated with the usual coordinate angles 6 and 4>- The 
helicity operator Jr = e.^ ■ J generates local rotations 
about the radial direction. In the helicity basis, compo- 
nents of 3-dimensional tensors have definite spin weight, 
which can be determined by counting the number of -|- 
and — index values. For example, consider a tensor com- 
ponent = T = e+ • r • e+. It is a spin-weight-2 
function and can be expanded into the complete set of 
spin-weight-2 spherical harmonics 2^m- Another exam- 
ple is the tensor component T_,. = e_ ■ T • e^, which has 
a spin weight of —1. In the present context, the elec- 
tric and magnetic fields have spin-weight-0 components 
Er = e-r ■ E, Br = Br ■ B and components E± = e± ■ E, 
B± = e± ■ B, which have a spin weight of ±1. 

In general, harmonics of a given weight satisfy comple- 
tion and orthogonality relations, 

^ sY*, (fi) sYjm m = s{n-n'), (Al) 

jm 

J sY*M ,Yj,m' {^) dn = 8jj,5mm' ■ (A2) 

More generally, they obey 

V. Y- - / (2ji-H)(2i2"+l) ' 

Sl-IjimiS2-'j2TO2 — Y 47r(2j3-M) 

83^3 "13 

X (ili2(-Sl)(-S2)|i3(-S3)) 

X (Ji32mim2\hmz) s^Yj^ms , (A3) 

where {jij2'm\m2\j3,m^) are Clebsch-Gordan coefficients. 
Note that orthonormality and completion do not extend 
across harmonics of different spin weight. 

The spin-weighted harmonics can be generated for 

the familiar spin-weight-zero harmonics o^m = Yjm 
using covariant-angular-momentum operators Jo = 
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^ieabcX^'^'^ + Sa^ wherc V represent covariant deriva- 
tives and Sa are covariant spin operators. The spin op- 
erators act on tensors in a manner similar to a connec- 
tion. For example, operating on a tensor T^o) we get 
SaT\ = Sl^T\ - StT\, where 5^, = ie^b- The 
totally antisymmetric tensor has nonzero components 
s+r- = —£"'"'"" = i in the helicity basis. In general, the 
components J± = e±- J raise/lower the spin weight of a 
function. Specifically acting on the spherical harmonics, 
we get the ladder relations 

J± sYjm = -^^{j{j + l)-s{s±l))s±lYjm ■ (A4) 

Successive operations generate arbitrary sYjm in terms 



of oXy m- 

Finally, there is a useful relationship between the co- 
variant derivatives and the angular-momentum operators 
in the helicity basis: 

= d/dr , 

V± = ±( J± - 5±)/r . (A5) 

These identities help in the reduction of differential ten- 
sor equations into radial and angular parts. For ex- 
ample, the divergence of vector F° becomes VaV^" = 

(I? + f)^'- + ~ •^-^+)' usi'^g properties of 

the spin operator. 
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